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ABSTRACT 

The on-orbit radiometric response calibration of the VISible/Near InfraRed (VISNIR) and the Short-Wave InfraRed 
(SWIR) bands of the Visible/Infrared Imager/Radiometer Suite (VIIRS) aboard the Suomi National Polar-orbiting 
Partnership (NPP) satellite is carried out through a Solar Diffuser (SD). The transmittance of the SD screen and the 
SD’s Bidirectional Reflectance Distribution Function (BRDF) are measured before launch and tabulated, allowing 
the VIIRS sensor aperture spectral radiance to be accurately determined. The radiometric response of a detector is 
described by a quadratic polynomial of the detector’s digital number (dn). The coefficients were determined before 
launch. Once on orbit, the coefficients are assumed to change by a common factor: the F-factor. The radiance 
scattered from the SD allows the determination of the F-factor. In this Proceeding, we describe the methodology and 
the associated algorithms in the determination of the F-factors and discuss the results. 
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1. INTRODUCTION 

The Visible Infrared Imaging Radiometer Suite (VIIRS) aboard the Suomi NPP satellite (launched on October 
28, 201 1) is the next generation polar-orbiting remote sensing whisk-broom scanner to the MODerate Resolution 
Imaging Spectroradiometer (MODIS), providing imaging capabilities in 22 spectral bands over a wavelength range 
of 0.4 pm to 12.0 pm . The design specified and the measured band central wavelengths and widths can be found in 
Ref. 1. Although carefully calibrated before launch, the on-orbit VIIRS sensor’s radiometric response can change 
substantially with time as seen in MODIS 2 . Consequently, on-orbit radiometric calibration for the VIIRS sensor is 
critical for generating high quality weather and environmental products. In this Proceeding, we give the 
methodology and the results of the on-orbit calibration of the Reflective Solar Bands (RSB). Our results are used to 
validate the results from the Aerospace Corp. 

2. METHOD FOR THE ON-ORBIT RADIOMETRIC RESPONSE CALIBRATION 

The primary on-orbit radiometric calibration in the RSB for NPP VIIRS is through the onboard Solar Diffuser 
(SD). Fig. 1 is a schematic showing the onboard RSB calibration physical components. When fully illuminated by 
the Sun through a pinhole attenuation screen, the SD provides a known spectral radiance source over a wide 
spectrum covering the visible and near infrared wavelengths by scattering the incident solar radiation. As a result, 
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viewing the fully illuminated SD allows the determination of a detector’s response coefficient at the particular 
spectral radiance level provided by the SD. The validity of the coefficient is then extended to other radiance levels. 
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Fig. 1. A schematic shows the observation of the SD scattered solar radiation. 

We denote the response coefficient as F, used as a multiplication factor in the following equation 

L sd = Fx(c 0 + q *dn + c 2 *dn 2 ^ , (1) 

where dn is a detector’s digital number with background subtracted. In Eq. (1), (c 0 , q, c 2 ) are determined through 
linear interpolations in the temperature of the electronics and in the temperature of the VISNIR focal plane for the 
VISNIR bands and the Optical Mechanics Module for the SWIR bands using prelaunch determined (c 0 , q, c 2 ) 3 at 
three temperature plateaus. L SD is the solar spectral radiance scattered from the SD, computed by 


Lsd = °;' ar T al P ° Wer x% x cos d SD _ sun x BRDF xRVS . (2) 

4x*dniRS-sun 

Eq. (2) is critical to the computation of the F-factor defined in Eq. (1) and needs careful explanations. In Eq. (2), 
® solar spectral power * s the solar spectral power at a distance of one astronomical unit, a wavelength-dependent 

predetermined time-independent quantity (the time independence is an approximation). r SAS is the transmittance of 
the pinhole screen for the SD, to allow the SD to provide spectral radiances close to the typical radiances for the 
RSB bands. Due to the pinholes, r SAS is dependent on the incident solar angle and can also be weakly dependent on 
the wavelength due to pinhole inner surface total reflection. As shown in Fig. 2, the incident solar angle components 
are quantified by a projected vertical angle (j) v in the XZ plane and a horizontal angle (j) H in the XY plane in VIIRS 
coordinate system. The X axis is along the longitudinal direction of VIIRS and the Z axis points to the EV port. 
Mathematically, the angles are computed by 

(j) H = -arctan(y/x) , (3) 

and 
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(fry = arctan(z/x) . 


( 4 ) 



Fig. 2. Definitions of the projected vertical and horizontal angles. 

t sas is determined before launch. 0 SD _ sun is the angle between the normal of the SD surface and the SD-Sun 
vector. BRDF is the Bidirectional Reflectance Distribution Function of the SD, a function of both the incident solar 
angle and the VIIRS sensor’s SD view angle as well as wavelength and time. Since the VIIRS sensor’s SD view 
angle is time invariant, BRDF may be expressed as 

BRDF = h(/ l, t-t 0 )x BRDF(/ 1, t = t 0 ; (fr H , (fry ) , (5) 

where H{A, t-t 0 ) represents the deterioration of the SD due to solar radiation, t 0 is the starting time when the SD is 
initially exposed to the Sun, and BRDF{z l, t = t 0 ; (fr H , (fry) is the initial BRDF of the SD, determined before launch. 

0 ) is determined by an onboard optical device which views the Sun (through a pinhole screen) and the SD 
at almost the same time. RVS (Response Versus Scan) is the reflectivity of the Half- Angle-Mirror (HAM) of VIIRS, 
an incident angle dependent quantity 4 . The HAM is an integral part of the sensor, residing right after the Rotating 
Telescope Assembly (RTA) in the optical path. RVS is weakly wavelength dependent for the RSB. Note that unlike 
other components in the sensor’s optical path, the HAM’s reflectivity is treated independently because as the HAM 
scans the EV port, the incident angle varies. 

Before we can use Eq. (1) to compute the F-factor, we need to pay attention to a detector’s ability to convert 
incident photons to electrons over a wide range of wavelengths. This ability is quantified by the Relative Spectral 
Response function 1 , denoted by RSR(A, b) with B denoting a particular band. With RSR , Eq. (1) is modified to be 

p \RSR(X, B)x L SD {X)dX 

(c 0 + q * dn + c 2 * dn 2 jx J RSR(A, B)dA 

3. SWEET SPOT DETERMINATION 

When the SD is fully illuminated by the Sun, seen by the VIIRS detectors, the solar vector angle components 
(fr H and (fry are said to be in the sweet spot. VIIRS is such designed that when at the sweet spot, the Z axis points at 
part of the Earth which is dark, avoiding complicated earth shine through the EV port in the F-factor determination. 
In addition, the louvers of the SD screen are designed to reduce earth shine on the SD. The sweet spot is determined 
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by plotting the DN (not background subtracted dn) of the detectors versus the projected vertical solar angle (j) v . A 
set of plots of DN versus (j) v is shown in Figs. 3a-b for orbits 154 (2011-11-08) and 3272 (2012-06-15). 
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Fig. 3. Plots of DN(SD) (not background subtracted) versus (j) v for orbit 154 (a) and orbit 3272 (b) for the first and the last 

detectors of bands II and Ml respectively, with (f) H ^28.0 (orbit 154) and (j) H ^15.5 (orbit 3272), respectively. The lower 
curves for band Ml are for the low gain state and the higher curves for the high gain state. 

We can see in Figs. 3a-b that when (j) v is greater than 19 degrees, DN drops very quickly with increasing angle, 
a result of earth shadow covering the SD; and when (j) v is less than about 6 degrees, DN drops quickly with 
decreasing angle, due to the moving out-of-detector- sight of the solar illumination on the SD at smaller (j)y . From 
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Fig. 3, the sweet spot is chosen to be (j) v in the region of [14.0, 18.0] degrees (the lower limit of the sweet sport can 
be relaxed to be a smaller angle). The sweet spot is much less sensitive to the horizontal angle (j) H evidenced by the 
similarity of the plots although the corresponding horizontal angles are quite different (at orbit 154 (j) H « 28.0 and 
at orbit 3272 (/> H « 15.5 ) . The lower curves for band Ml in Figs. 3a-b are for the low gain state. 

4. RESULTS 

After the degassing phase, NPP VIIRS began operation mode on November 8, 201 1 on orbit 154, allowing the 
start of the F-factor determination for the Visible and Near InfraRed (VISNIR) bands in the RSB. The F-factor 
determination for the Short Wave InfraRed (SWIR) bands in the RSB started after the cryo-door opened a few days 
prior to orbit 1201 to cool down the SWIR band focal plane. 

For each orbit, there are a few dozen scans satisfying the sweet spot criterion, with each of the scans lasting 
nominally 1.7864 seconds. We compute the F-factor for each scan and then compute the simple average per orbit 
through a repeated 4-sigma outlier rejection procedure. We compute the F-factor for each detector in each of the 
fourteen RSB bands and at both HAM sides and in both high and low gains for dual gain bands. To give a few 
examples of the computed F-factors, in Fig. 4, we plot per orbit simple-averaged F-factors versus orbit number (each 
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Fig. 4. F-factors versus orbit number for the mid-detectors of RSB bands at HAM side A and high gain for dual gain bands. 
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orbit lasts about 101 minutes) for the mid-detector (detector 8 for the M bands and 16 for the I bands) of each of the 
RSB bands at HAM side A and high gain for dual gain bands. A major feature in Fig. 4 is that the F-factors for the 
Ml band are about 1.3, much larger than one. The reason is due to that in the prelaunch calibration to determine 
(c 0 , c x , c 2 ) , the RSR in the out-of-band region was not considered. The unusually large out-of-band RSR at 
wavelengths larger than that for the central peak for the Ml band, coupled with a radiance source redder than the 
Sun, leads to a computed spectral radiance smaller than what it should be. Consequently, c x for the Ml band is 
smaller than what it should be, yielding a large F-factor. The other major feature in Fig. 4 is that for bands II, 12, 

M5, M6, and M7, the computed F-factors increase quickly with increasing orbit numbers after orbit 344 during 
which the NADIR door was opened, indicating degradation of optical components along the sensor’s optical path. 
The F-factors for the SWIR bands also show such increase. Speculation on the cause of the degradation includes 
solar radiation damage. The radiation is mainly from the reflection of the solar radiation by the Earth through the EV 
port and in a small portion from the scattered solar radiation by the SD. 

To help identify whether solar radiation is the root cause, the RTA was stowed for a few times, pointing away 
from the EV port. During the stows, no solar radiation entered the RTA. We plot in Fig. 5a, the VISNIR band mid- 
detector per orbit averaged F-factors normalized to the respective values at orbit 154 versus equivalent EV port solar 
radiation exposure time. In Fig. 5b, we plot the SWIR band mid-detector per orbit averaged F-factors normalized to 
the respective values at orbit 1201 versus the exposure time. The plots show smooth curves, revealing that solar 
radiation is very likely the cause of the degradation. 



Fig. 5. Per orbit averaged normalized F(mid det)/F(mid det; orbit 154) versus EV port solar radiation exposure time for the (a) 

VISNIR and (b) SWIR bands 
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Reflectance and surface contamination measurements on the Telescope Witness Mirror (TWM) by the 
Aerospace Corp. indicate that the Rotating Telescope Assembly may be contaminated with tungsten oxides on its 
four reflective mirror surfaces. It is known that tungsten oxides, after exposure to UV radiation, become light 
absorbent over NIR and SWIR wavelengths 5 . Mathematical models to account for the observed F-factors, based on 
RTA mirror surface contamination, have been developed. The models and the applications of the models to address 
issues such as projected and smoothed F-factors, projected signal-to-noise ratio, and the impact of the RTA mirror 
darkening effect on the RSR will be described elsewhere. 

In addition to the unusually large F-factors for band Ml and the quick increase of the F-factors with increasing 
orbit numbers after orbit 344 for bands II, 12, M5, M6, and M7, there are a few noticeable features for the F-factors. 
To help visualize the features, we plot in Figs. 6a-b, per orbit averaged normalized F-factors versus orbit number. 
First, the large deviations from the trend around orbits 1563 to 1577 are due to the yaw maneuvers of the spacecraft. 
In addition, for the SWIR bands, the discontinuities around orbits 2101 to 2133 and the unusually large F-factors 
right after orbit 2133 are due to the Spacecraft Control computer reset. The reason why the reset caused the behavior 
of the F-factors is under investigation. We also notice that many of the peaks and valleys in the F-factors are in sync 
across the bands. In particular, a hump is observed over orbits 1800 to 1900 for all bands, likely due to a solar storm. 



Fig. 6. Per orbit averaged normalized F(mid det)/F(mid det; orbit 154) versus orbit number for the (a) VISNIR and (b) SWIR 

bands 

Finally, as revealed by Eqs. (2) and (6), there are many sources of error in the computed F-factor. For example, 
an error in t-t 0 ) would propagate into the computed F-factor. Quantitative descriptions of the uncertainty of 
the computed F-factor will be addressed elsewhere. 
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5. SUMMARY 


In this Proceeding, we described the methodology for the on-orbit radiometric response calibration through the 
SD for the Suomi NPP VIIRS RSB bands. We computed the F-factors for each detector in the RSB bands at both 
HAM sides and in both high and low gains for dual gain bands. The F-factors for the NIR and the SWIR bands 
increase quickly with time right after the NADIR door was opened in orbit 344, demonstrating unexpected 
degradation of some of the components along the sensor’s optical path. Plots of the F-factors versus the EV port 
solar radiation exposure time yield smooth curves, indicating that the source of the degradation is likely the solar 
radiation, consistent with the contamination of tungsten oxides on the RTA reflective mirror surfaces. The computed 
F-factors were used to validate the results from the Aerospace Corp. 
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